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PARKING FUNCTIONS AND TREE INVERSIONS REVISITED 


PETAR GAYDAROV AND SAM HOPKINS 


Abstract. Kreweras proved that the reversed sum enumerator for parking func¬ 
tions of length n is equal to the inversion enumerator for labeled trees on n + 1 
vertices. Recently, Perkinson, Yang, and Yu gave a bijective proof of this equality 
that moreover generalizes to graphical parking functions. Using a depth-first search 
variant of Dhar’s burning algorithm they proved that the reversed sum enumerator 
for G-parking functions equals the K-number enumerator for spanning trees of G. 
The K-number is a kind of generalized tree inversion number originally defined by 
Gessel. We extend the work of Perkinson-Yang-Yu to what are referred to as “gener¬ 
alized parking functions” in the literature, but which we prefer to call vector parking 
functions because they depend on a choice of vector x £ N n . Specifically, we give an 
expression for the reversed sum enumerator for x-parking functions in terms of in¬ 
versions in rooted plane trees with respect to certain admissible vertex orders. Along 
the way we clarify the relationship between graphical and vector parking functions. 


1. Introduction 

Set N := {0,1,2,...}. Parking functions are certain sequences in N n that originally 
arose in the study of hashing functions [23] , Their name comes from the following 
description in terms of parking preferences for cars: 

Imagine that n parking spots labeled 0 through n — 1 are arranged in 
order on a linear street. Cars C\,... ,C n approach the spaces in order. 

Car Ci prefers spot aq, which means that he will drive until he reaches 
that spot and park there it if it is unoccupied. If his preferred spot 
is occupied he will continue driving until he comes to an unoccupied 
spot and park in this unoccupied spot. If his preferred spot and all the 
spots after it are occupied, then he cannot park. A parking function of 
length n is sequence (or,... , a n ) £ N n of parking preferences so that 
all the cars can park. 

An equivalent, less whimsical, definition of parking function is the following. 

Definition 1.1. A parking function of length n is a sequence (oq,. .., a n ) £ N n whose 
weakly increasing rearrangement ccq < cu 2 < ■ ■ ■ < aq n satisfies oq. < j — 1 for 
all 1 < j < n. We denote the set of parking functions of length n by PF(n). 

Parking functions have appeared in various of areas of pure mathematics. We will not 
attempt any kind of survey here except to say that, beyond their inherent combinatorial 
interest, parking functions arise in the study of hyperplane arrangements 001 Lecture 6] 
and the theory of diagonal harmonics PH Chapter 5]. For more background on parking 
functions in general, see the recent survey by Yan [45] . 
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Parking functions are intimately related to labeled trees. For instance, it is well- 
known that the number of parking functions of length n is (n + l) n_1 [241137] . Cayley’s 
formula tells us that this is also the number of labeled trees onn + 1 vertices. Indeed, 
there are many explicit bijections between parking functions and labeled trees [581 [T3] . 
But more is true. The sum of a parking function a = (cci,... ,a n ) £ PF(n) is defined 
to be sum(a) := Y^i=\ a i an d the reversed sum is rsum(a) := ( 2 ) — sum(a). The 
reversed sum is defined so that the parking functions of maximal sum have reversed 
sum 0. The reversed sum is also called the “area of a major sequence” in [4] and 
“total displacement” in [23] but we follow the terminology used in [28] . Let us call 
the expression X^aGPF(n) q ,rsum( ' C d the reversed sum enumerator for parking functions 
of length n. Kreweras [ 26] established that the reversed sum enumerator for parking 
functions of length n is equal to the inversion enumerator for labeled trees on n + 1 
vertices. That is, Kreweras proved that 

(1.1) g rsum ( Q ) = ^ g inv ( T ) 

aePF(n) T 

where the sum over all labeled trees T on n + 1 vertices and inv(T) is the inversion 
number of T, a certain natural statistic that generalizes the inversion number of a 
permutation. Specifically, for T a labeled tree with vertices 0,1,2,... , n (which we will 
always consider rooted at 0) an inversion of T is a pair (i, j) with 1 < i < j < n such 
that the unique path from i to 0 passes through j. The inversion number inv(T) is the 
number of inversions of T. 

on 3 vertices and their inversion 


o 0 
u 2 
" 1 

inv(T) = 1 

On the other hand, the three parking functions of length 2 are (0, 0), (1, 0), and (0,1), 
which have reversed sums of 1, 0, and 0 respectively. So m holds in this case. 

For more on the inversion enumerator for labeled trees see m and [18]. Kreweras’s 
proof of (11.11) was via a recursively defined bijection between parking functions and 
trees that preserves the requisite statistics. Stanley stated the problem of finding a 
nonrecursive bijective proof of (11.11) as [4D[ Lecture 6, Exercise 4]. Various authors [59l 
[20] devised such a nonrecursive bijection. Recently, Perkinson, Yang, and Yu [35] found 
an elegant bijective proof of (11.11) that moreover generalizes naturally to the setting of 
graphical parking functions. To keep them distinct from their generalizations, let us 
refer to the regular parking functions defined above as classical parking functions. 

Graphical parking functions were first given that name by Postnikov and Shaprio [36] 
but have been studied for longer in the context of the abelian sandpile model where 
they are essentially the same as what are called superstable configurations (see m 
Definition 4.3] or [ 34l §2.3]). The superstable configurations are dual to the recurrent 


Example 1.2. Set n := 2. The three labeled trees 
numbers are the following: 


0 

«4*Y 

1 2 
inv(T) = 0 


0 
1 
2 

inv(T) = 0 
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configurations (see (22, Definition 2.11] or [341 Definition 2.10]; these recurrent configu¬ 
rations are also sometimes called critical configurations as in ED- The set of recurrent 
configurations is a more natural object to study from the perspective of statistical me¬ 
chanics, the field out of which the abelian sandpile model originally arose. Indeed, the 
abelian sandpile model was first described in the pioneering work of Bak, Tang, and 
Wiesenfeld |2j, who were interested in modeling self-organized criticality in physical 
systems. Dhar [12] subsequently put the model into a more general graphical frame¬ 
work and developed much of the basic theory. Let us also note that graphical parking 
functions appear under the name of vq- reduced divisors in Riemann-Roch theory for 
graphs [3j. As such they are related to a deep analogy between finite graphs and Rie- 
rnann surfaces. For more background on sandpile theory, consult the short survey [29] 
or the upcoming book Ell- 

Graphical parking functions depend on a choice of graph G (as well as a choice of 
root vertex that we will ignore in this introduction). We denote the set of G-parking 
functions by PF(G). For a € PF(G) we have similar notions of sum sum(a) (which in 
this context is often called degree ) and reversed sum rsum(a). G-parking functions serve 
as coset representatives for the cokernel of the reduced Laplacian A(G) of the graph G. 
Therefore, by Kirchhoff’s Matrix-Tree Theorem, the number of G-parking functions is 
equal to the number of spanning trees of G. It is thus natural to ask for an extension 
of Kreweras’s result (11.11) to this graphical setting which says that the reversed sum 
enumerator for G-parking functions is some kind of “inversion enumerator” for spanning 
trees of G. 

Perkinson-Yang-Yu [35 j provided such an extension. Specifically, they gave a bijec- 
tive proof that 

(1.2) ^ grsum(q) _ ^ g«(G,T) 

aePF(G) T 

for any simple graph G, where the sum is over all spanning trees T of G and k(G, T) is 
a kind of generalized inversion number originally defined by Gessel m- Importantly, 
when G = K n+ \ is the complete graph, equation (11.21) recovers (11.11) . The proof of (11.21) 
in [35] is based on a depth-first search (DFS) variant of Dhar’s burning algorithm ~[~2i 
Eg. The idea of using depth-first search to study tree inversions goes back at least to 
the fundamental work of Gessel and Wang m • Depth-first search appeared again in 
the paper m in which Gessel defined the K-number statistic. And depth-first search, 
as well as the related neighbor-first search, were also used later by Gessel and Sagan HU 
to relate certain spanning tree activities (including generalized inversion numbers) to 
the Tutte polynomial of a graph and to classical parking functions. The main novelty 
of [35] was the combination of depth-first search with Dhar’s burning algorithm. (But 
note that similar ideas appeared in [9, §3] and especially in [ 251 §4], which however used 
neighbor-first search rather than depth-first search per se.) Dhar’s burning algorithm 
is one of the most important tools in the study of the abelian sandpile model: it is a 
linear time algorithm for deciding whether a given sequence is a G-parking function 
that in some sense generalizes the “parking procedure” mentioned at the beginning of 
this section. Another variant of the burning algorithm, due to Cori and Le Borgne m, 
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proves that the reversed sum enumerator for G-parking functions is also the generating 
function for spanning trees of G by external activity (a different tree statistic that 
depends on a total edge order rather than a total vertex order). 

Remark 1.3. There is a natural embedding of the symmetric group & n on n letters 
into the set of labeled trees on n + 1 vertices given by a eA T a where T a is 

0 o-(l) a (2) cr(n) 

a a-«-- 9 

Observe that inv(T CT ) = inv(cr), the usual permutation inversion number. Pushing T a 
through the bijection of Perkinson-Yang-Yu [35] yields not the most obvious way of 
representing a in PF(n) (“one-line notation”) but rather a version of the code of a 
(see [3TJ §1.3]). Indeed, the identity (II.ip should be seen as a generalization of the very 
classical result 

[»]„! = E “‘ nrM 

cr£&n 

where [n\ q \ := [n] q ■ [n - 1] 9 • • • [2], ■ [1], and [k] q := = 1 + q + q 2 -\ -h q fe_1 are 

the standard q-factorials and q-numbers. This offers some impression of the subtlety 
of the DFS-burning bijection. 


There is another generalization of classical parking function that goes under the 
unfortunately vague name of generalized parking function [33] [37]. Whereas graphical 
parking functions depend on a choice of graph G , generalized parking functions depend 
on a choice of vector x = (x \,... ,x n ) £ N n of nonnegative integers. As such we will 
call them x- parking functions or vector parking functions. 


Definition 1.4. Let x = (x\,... ,x n ) £ N n be some vector of nonnegative integers. 
An x- parking function is a sequence (aq, ...,a n ) € N n whose weakly increasing re¬ 
arrangement ct^ < oti 2 < ■ ■ ■ < a.i n satisfies < (Yjk =i x k) ~ 1 for all 1 < j < n. We 
denote the set of x-parking functions by PF(x). 


Observe that PF(n) = PF(x) for x = (1,1,... , 1). Also note that rational parking 
functions, which because of their connections to algebra and geometry (especially the 
representation theory of rational Cherednik algebras 0) have been much studied as of 
late [mm], are a special case of vector parking functions. But vector parking functions 
are less structured objects than rational parking functions: for instance, it is easy to 
count rational parking functions (see [1] Corollary 4]) but apparently not so easy to 
count arbitrary vector parking functions. The best general formula for enumerating x- 
parking functions is due to Pitman and Stanley m Theorems 1 and 11] and appeared 
in the paper in which they introduced x-parking functions. In order to state their result 
we need a little more notation: for all n £ N, set 


r(n) 


j 

{(7i,---,7n) €N n : E 7 i > j for all 1 < i < n 

i= 1 


n 

1 and ^2 yj = n}. 
2—1 
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It is well-known that ffT(ri) = C n where C n := ■^qrx( 2 f ") is the nth Catalan number; 
see e.g. @2 §2 Exercise 86]. 


Theorem 1.5 (Pitman-Stanley). For any x = (xi,...,x n ) G N n , the number of x- 
parking functions is the following polynomial P n (x i,..., x n ) in the variables xi ,..., x n : 


Pn(x 1j■ ■ ■ j Xn) 


E 


^'Qil+l^Q'2 + 1 ’ * * ^OLn -\-1 


(ai,...,a n )GPF(n) 


E 

(7l.---,7n)er(n) 


n\ 

7 i! 72 !• • 


In 


™7l ™72 
| X 1 x 2 


rpln 


Vector parking functions have been studied in many subsequent papers, especially 
by Yan and her collaborators. In [3?] Yan gave nice enumerative formulas for several 
special classes of x-parking functions. In [44], building on earlier work [46], Yan related 
the reversed sum enumerator for a special class of x-parking functions to tree inversions. 
In [SI] an d m Kung and Yan studied the connection between x-parking functions 
and Goncarov polynomials. Via this connection they offered linear recursions, a shift 
formula, and a perturbation formula for the number of x-parking functions, as well as 
a formula for the expected sum of a random x-parking function. They also investigated 
the reversed sum enumerator for PF(x). For x = (xi,... ,x n ) G N n , define the sum 
of a vector parking function a = (aq,... ,a n ) € PF(x) to be sum(a) := Y17= l ai an< ^ 
the reversed sum to be rsum(a) := (X^=i( n + 1 — i)xi) ~n — sum(a). Again, reversed 
sum is defined so that x-parking functions of maximal sum have reversed sum 0. Kung 
and Yan [28, Theorem 6.2] established the following formula for the reversed sum 
enumerator for PF(x) (and they also explained how it follows from the arguments of 
Pitman-Stanley ji3|). 


Theorem 1.6 (Pitman-Stanley and Kung-Yan). For any x = (xi,... ,x n ) G N n , the 
reversed sum enumerator for x -parking functions is 


£ ^ rSUm(a) = q^ {n+1 - i)xi) ~ n Pn([xi }i/ q ,q- xi [x 2 ]i/ q , ■ ■ .,Y I1 - I2 -- I "- 1 [xn]i/,) 

agPF(x) 


- E 

(71,-,7n)er(n) 


n\ 

7 i! 72 !•• 


_gE"=i 1 (7i+72-l—^ -+7i-i)xi+i 

7n! 



Since the best general enumerative formula for ^PF(x) (Theorem 11.511 involves a 
sum of Catalan many terms we cannot hope for a better formula for the reversed 
sum enumerator for PF(x) than one which also sums Catalan many terms. Thus in 
some sense Theorem 11.61 is totally satisfactory. But on the other hand, in light of the 
connections between tree inversions and the reversed sum enumerators for classical and 
graphical parking functions, we might hope that we could also express the reversed 
sum enumerator for vector parking functions in terms of tree inversions. In particular 
we could ask for an expression for the reversed sum enumerator for x-parking functions 
that evidently recaptures (11.11) when x := (1,1,..., 1) and Theorem 11.51 when q := 1. 
This is the problem we take up in the current paper. 
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We adapt the technique of Perkinson-Yang-Yu to give an expression for the reversed 
sum enumerator for x-parking functions in terms of tree inversions. In the special 
case where x = (o, b, b,... , b) £ N n we recover a result of Yan [33] that gives the 
reversed sum enumerator for PF(x) in terms of inversions in spanning trees of a certain 
multigraph (see Remark 13.41) . In this case the x-parking functions are exactly 

the graphical parking functions of K^+i- In fact, the first step of our proof is to 
extend in a straightforward way the DFS-burning algorithm of [35] to multigraphs (see 
Theorem 13.11) . But our result is more general than that of [44j because it applies to 
arbitrary x. In general vector parking functions are not graphical parking functions of 
any multigraph and vice-versa. We completely characterize the (small) overlap between 
graphical and vector parking functions in Theorem 12.51 

Our main result expresses the reversed sum enumerator for PF(x) as a sum over 
rooted plane trees. In order to state this result, let us briefly review notation for trees. 
A rooted tree is a tree together with a choice of distinguished root vertex. A rooted 
tree T comes with a partial order <t on its vertex set V (T) whereby i <t j if and only 
if the unique path in T from j to the root passes through i. If i <t j and i ^ j then 
we say i is an ancestor of j and j is a descendant of i. We use other familial language 
in the obvious way: so we say that j is a child of i, and i is the parent of j, if j is a 

descendant of i and no descendant of i is an ancestor of j. For a non-root vertex i, 

we use par T (j) to denote the (unique) parent of j. We say that i and j are siblings 
if par T ({) = par r (j). A rooted plane tree is an rooted tree drawn in the plane with 
parents drawn above children and with the children of each vertex linearly ordered from 
left to right. One fundamental combinatorial interpretation of the Catalan number C n 
is as the number of rooted plane trees on n + 1 vertices (see @21 Theorem 1.5.1 (iii)]). 
By giving the tree a depth-first labeling (or preorder as in @21 p. 10]) of its vertices, 
we have the following equivalent definition of rooted plane tree. 

Definition 1.7. For n > 1, a rooted plane tree onn + 1 vertices is a rooted tree T with 
vertex set V (T) := {0,1,..., n} and root 0 satisfying the following conditions: 

• i <r j implies i < j for all 1 < i, j < n; 

• i <t k implies i <t j for all 1 < i < j < k < n. 

We denote the set of rooted plane trees on n + 1 vertices by RPT(n + 1). 

To state our main result we need also to discuss vertex orders for rooted trees. For T a 
rooted tree and -< any total order on its non-root vertices, an inversion ofT with respect 
to Y is a pair (i. j) of non-root vertices with i <t j but j -< i. For T £ RPT(n-|- 1), an 
admissible vertex order of T is a total order -< on V ( T) — {0} for which 1 < i < j < n 
with i and j siblings implies j -< i. We denote the set of admissible vertex orders 
of T by AVO(T). Observe that the number of pairs (T, -<) with T £ RPT(n + 1) 
and -<£ AVO(T) is (n + l)”^ 1 , the number of labeled trees onn + 1 vertices. Indeed, 
given any labeled tree T' with vertices 0,1,... ,n there is a canonical way to associate 
such a pair (T, -<): perform a depth-first search of T' starting at 0 and always preferring 
to visit the vertex with the greatest label when presented with a choice; then, if a £ 6 n 
is the unique permutation recording the order 0, cr(l), <r(2),... , <r(n) in which we visited 
the vertices of T' in the course of this search, relabel T' by i i —> o'~ 1 (i) to obtain a rooted 


PARKING FUNCTIONS AND TREE INVERSIONS REVISITED 


7 


0 



V 




0 



Figure 1. An example of the bijective correspondence between labeled 
trees and rooted plane trees together with an admissible vertex order. 
On the left we have an arbitrary labeled tree T' on 9 vertices. On the 
right we have the corresponding T G RPT(9) and -<!G AVO(T). 


plane tree T; lastly, define -< by * -< j if and only if a{i) < a(j). See Figured] for an 
example of this bijective correspondence. Note crucially that under this correspondence 
inversions of T' correspond to inversions of T with respect to -<. 

The main result of this paper is a bijective proof of the following theorem. 


Theorem 1.8. For x = (aq,..., x n ) G N n we have 


E » 

o;GPF(x) 


rsum(a) _ 


E 


f 


rr parj'(i) + l 

E 9 ^ 


TeRPT(n+l) \^gAVO(T) 



outdeg T ( 2 —1) 


where outdeg T (i) := ff{j € V(T): i = par T (j)} is the number of children of i in T. 


Theorem 1 1. 81 does relate the reversed sum enumerator for x-parking functions to tree 
inversions, albeit in a slightly complicated way. (The dependence of the expression on 
the parent of i for an inversion (i, j) is a phenomenon going back to Gessel’s original 
definition [16] of his r- number.) It is worth comparing Theorems ll.6l and ll.81 They both 
express the reversed sum enumerator for x-parking functions as a sum of Catalan many 
“terms.” The “main parts” (products of (/-numbers) of these terms are the same in the 
two formulas; however, the “coefficients” in front of the main parts are different. Thus 
Theorems dn and 11.81 together imply some surprising equalities between expressions 
involving q and the variables x \,..., x n . Indeed, we can see these surprising equalities 
emerge already in the smallest nontrivial case of n := 2 , as we detail in the following 
example. 

Example 1.9. Let x := (xi,X 2 ) G N 2 . First let us compute the reversed sum enu¬ 
merator for x-parking functions according to the formula in Theorem 11.61 In this 

case T(2) = {(2,0), (1,1)}. For 7 = (2,0) we have Ya= 1(71 d-h 7 i~ i)xi+i = x 2 

and for 7 = (1,1) we have X)i=i(7i + ' —h 7 i — i)xi +1 = 0. Thus from Theorem 11.61 we 
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conclude that 

g r sum(a) = < f>[ Xl )2 q + 2[x 1 ] q [x 2 \ q 

aGPF(x) 


Now let us compute the reversed sum enumerator for x-parking functions according to 
Theorem 11.81 So we fill out the following table: 


T E RPT(3) 

-<E AVO(T) 

Inversions of T 
w.r.t. -< 

• r par T (i)+l 

j-<i 


0 





2 -« 1 

0 

0 

1 

2 




< 

0 




< 

> 1 

1 -< 2 

0 

0 

i 

2 





0 




< 

» 1 

2 ^ 1 

{( 1 , 2 )} 

X\ 

( 

2 





Thus from Theorem 11.81 we conclude that 

E 9”“ l “ ) = K + (l + 9* , )N,W,. 

oGPF(x) 

In fact wehave g X 2 [xi]|+2[xi] ? [x 2 ] ? = [xi ] 2 + (l+g xi )[xi] 9 [x 2 ] g for all x = (xi,x 2 ) E N 2 , 
so Theorems 11.61 and 11.81 agree in this case. 

Example 1.10. For n := 3, Theorems 11.61 and 11.81 together imply that 
q 2x2+X3 [x i]J + Zq X2+X3 [xi] 2 q [x 2 ] q + &f a [xi] 2 q [x 3 ] q + 3q X3 [xi] q [x 2 ] 2 q + 6HJ* 2 ] g [z 3 ] g = 
[xi]q + (1 + 2 Q X1 )[xi} 2 [x 2 ] q + (2 + q xi ) [sq] 2 [x 3 ] q + (1 + q xi + q 2xi )[x l\ q [x 2 ] 2 q 

+ (1 + q Xl + q X2 + q 2xi + q Xl+X2 + q 2xi+X2 )[xi] g [x 2 }q[x 3 }q 

for all x = (xi, x 2 , x 3 ) E N 3 , an equality which is true but not obvious. 

By substituting q := 1 in Theorem 11.81 we obtain the following corollary. 


Corollary 1.11. For all x = (xi,... ,x n ) E N n ; the number of x-parking fimctions is 


n\ 

TeRPT(n+l) 



outdeg T (i-l) \ 

-3 - . 

outdeg T (i — 1 )! I 


Corollary 11.111 and Theorem 11.51 are easily seen to be equivalent, as we explain in 
Remark 14.31 below. Moreover, the proof of Theorem II .51 appearing in [J3] is much more 
straightforward than the DFS-burning algorithm we use to establish Corollary 11.111 
We stress that the expression for the reversed sum enumerator in Theorem 11.81 is our 
main result. In particular we believe, as suggested by Examples 11.91 and 11.101 that 
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Theorems 11.61 and 11.81 are not equivalent in any obvious way and thus Theorem 1 1.8 1 is a 
genuinely new expression for the reversed sum enumerator, one which incorporates tree 
inversions. A feature of Theorem 11.81 which Theorem 11.61 lacks is that, thanks to the 
correspondence detailed above between labeled trees T' and pairs (T, -<) of rooted plane 
trees with admissible vertex orders, Theorem II.81 recovers the foundational result (11.11) 
of Kreweras when x := (1,1,..., 1). 

The moral of this paper is that while G-parking functions and x-parking functions 
are two largely orthogonal generalizations of classical parking functions, techniques 
from the former can fruitfully be applied to the latter, especially when it comes to the 
relationship between trees and parking functions. 

Acknowledgements: This research was carried out at MIT as part of the RSI summer 
mathematics research program for high-school students. The second author was the 
mentor of the first author. We thank David Perkinson, Alex Postnikov, and Richard 
Stanley for some helpful comments. We also thank an anonymous referee for very 
many helpful comments, especially with regards to the exposition of the paper and 
with pointers to the (vast) parking function literature. The second author was partially 
supported by NSF grant #1122374. 

2. The overlap between graphical and vector parking functions 

In this section we review graphical parking functions with the aim of understanding 
their overlap with vector parking functions. Beyond Definition 12.11 none of this section 
is necessary for the proof of our main result ('Theorem 1 1.81) . But Theorem l2.5l below does 
put our work in context because it shows that graphical and vector parking functions 
are genuinely different objects. 

Throughout this paper we will work with (finite) multigraphs , i.e., finite, undirected 
graphs for which multiple edges, but not loops, are allowed. Formally, we view a 
multigraph G as consisting of its vertex set V (G) together with its edge-weight func¬ 
tion log '■ However, we will also often think about G as consisting of its 

vertex set V(G) together with its multiset E(G) of edges: so E(G) contains wc(e) 
many copies of e for each e £ (^^). By ffE(G) := X^ ee (V(G)\ wc(e) we will mean the 

number of edges of G counted with multiplicity. Any multigraph G under consideration 
will come with a distinguished root vertex. All constructions that follow will implicitly 
depend on this choice of root but we will suppress the dependence in our notation. We 
use V ( G ) to denote the set of non-root vertices of G. For the rest of this section, G is 
a fixed multigraph with V{G) := {0,1,... ,n} (where n > 1) and root 0. 

Definition 2.1. A G-parking function is a formal sum a = JT G y (G) «i[*]eNR(G) of 

non-root vertices such that for any 0 # U C V(G) there is i £ U with a.{ < deg®(z) — 1 
where deg^(i) := Y,jev(G)-u j})- We identify the formal sum a = J2 iG v{G) 

with the sequence a = (aq,..., a n ) £ N n when V(G) = {1,..., n}. We denote the set 
of G-parking functions by PF(G). The sum of a £ PF(G) is sum(a) := Ya =i and 
the reversed sum is rsum(a) := ffE(G) — n — sum(a). 
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We say that G is connected if for each 0 / U C V ( G ) there is i G U with deg G (z) > 1. 
Note that PF(G) ^ 0 if and only if G is connected. Also note that PF(G) = PF(n) 
when G = K n +1 is the complete graph, i.e., when wc(e) = 1 for all e G So 

indeed graphical parking functions are a generalization of classical parking functions. 
Let us now recall some other basic facts about graphical parking functions. The first 
fact, which we also mentioned in Section |T] says that the number of G-parking functions 
is the number of spanning trees of G (weighted by edge multiplicities). A tree T is a 
connected multigraph with ffE(T) = ffV{T) — 1. We will always consider a tree 
to be rooted at 0 £ V(T). Therefore every tree T comes with a partial order <t 
on V (T) as described in Section [I] A spanning subgraph H of G is a multigraph H 
with V(H) = V(G) and E{H) C E{G). A spanning tree T of G is a spanning subgraph 
of G which is a tree. We use SPT(G) to denote the set of spanning trees of G. 

Fact 2.2. We have 

#PF(G)= ]T n u G (e). 

TeSPT(G) eGE(T) 

Thus via the Matrix-Tree Theorem (see |32l §5]J we conclude that #PF(G) is also equal 
to the determinant of the reduced Laplacian A(G) of G. 

Proof. This is very well-known; see for instance m Theorem 2.1]. At any rate, it will 
follow from Theorem o below by substituting q := 1. □ 

The next fact concerns maximal G-parking functions. The monoid N n has a natural 
partial order whereby for a = (aq ,... ,a n ), a' = (a \...., a' n ) G N n we write a < a' if 
and only if at < ol i for all 1 < i < n. A maximal G-parking function is a G-parking 
function that is maximal among G-parking functions with respect to this partial order. 
We use MPF(G) for the set of maximal G-parking functions. 

Fact 2.3. Fora G N n , we have a G PF(G) if and only if a < a' for some a' G MPF(G). 

Proof. This is an easy consequence of Definition 12.11 □ 

The final fact concerns another, extremely useful, combinatorial model for MPF(G) 
in terms of orientations. An orientation of G is a subset O C V(G ) 2 such that 

• (■ i , i) £ O for any z; 

• (i,j) G O implies w G ({i,j}) > 1 for all {i,j} G ( l ( 2 G) ); 

• for each {i,j} G ( V( 2 G) ) with uj G ({i,j}) > 1, #{(z, j), ( j,i )} (10 = 1. 

Given an orientation O of G and i G V(G), the indegree of i with respect to O 
is indeg c ,(i) := Yla j)eo w g({l j})- We say i G V(G) is a source of O if indeg c ,(i) = 0 
and we say it is a sink of O if indeg 0 (z) = deg G (z). Here the degree of i in G 
is deg G (z) := Yli^jev(G) w G({b j}) = c l e g{j}(0- Finally, we say O is acyclic if there does 
not exist a sequence (zi,^ 2 ), (* 2 ^ 3 ), • • •, (jkPk+i) with (ij,ij+ 1 ) G O for all 1 < j < k 
and i\ = ik+i- Let A(G) be the set of acyclic orientations of G with unique source 0. 

Fact 2.4. The map O i-A (indeg c ,(l) — 1, indeg G (2) — 1,... , indeg c ,(n) — 1) is a bijection 
between A{G) and MPF(G). 
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Proof. This is also well-known among those who study graphical parking functions; 
see Theorem 3.1]. The proof uses Dhar’s burning algorithm. □ 


The symmetric group & n acts on N n by cr (a) := (a^-icn,... for a G 6 n 

and a = (or, . .. ,a n ) € N n . We will now classify the multigraphs G for which PF(G) 
is invariant under this action of 6 n ; an immediate consequence will be a classification 
of those multigraphs whose parking functions are also the parking functions of some 
vector. To that end let us define some families of multigraphs. A cycle C is a connected 
multigraph with degc<(i) = 2 for all i G V(C). For a > 1 let us say G is an a-cycle 
if there is some cycle C with V{C) = V(G) and ojc(e) = a ■ u>c(e ) for all eG (' ^). 
Similarly, for a > 1 let us say G is an a-tree if there is some tree T with V(T) = V(G) 
and wc(e) = a-ujT(e) for all e G (' If'- 1 ). Finally, for a, b > 1 let \ be the multigraph 
with vertex set V(K%\) := { 0 , 1 ,... ,n} and edge-weight function 


({*, j}) 


a if i = 0 or j = 0 
b otherwise. 


Theorem 2.5. IfPF(G) is invariant under the action of <3 n on N n then one of the 
following holds: 

n —1 

• G is an a-tree and PF(G) = PF((a, 0,0,..., 0)) for some a > 1; 

n—2 

/- /V -s 

• G is an a-cycle and PF(G) = PF((o, 0, 0,... , 0, a)) for some a > 1; 

71—1 

• G = K^i and PF(G) = PF((a, b, b, ... , b )) for some a, b > 1 . 

On the other hand, ifPF(G) is not invariant under the action of the symmetric group 
then PF(G) ^ PF(x) for any xGff". 

Proof. From Definition [LT] it immediately follows that for all x G N n , PF(x) is invariant 
under the action of the symmetric group so the last sentence is certainly true. 

So let us consider when PF(G) can be invariant under & n . We claim that PF(G) is 
invariant under the action of the symmetric group if and only if MPF(G) is invariant 
under the action of the symmetric group. Fact 12.41 implies that all maximal G-parking 
functions have the same sum, which means that if a G MPF(G) and a (a) G PF(G) 
then a (a) G MPF(G). The other direction follows from Fact 12.31 

From now on we focus on maximal G-parking functions. We want to show that 
if MPF(G) is invariant under the action of the symmetric group then G is an a-tree 
or G is an a-cycle or G = We prove this by induction on n := ffV(G) — 1. 

The base case of n = 1 is trivial. So suppose that n > 2 and MPF(G) is invariant 
under & n . For i G V(G) let G — i denote the multigraph obtained by deleting i: 
this is the multigraph on vertex set V(G) — {i} and with edge-weight function equal 
to ujg restricted to ( v . A cut vertex of G is a vertex for which G — i is not 

connected. Define 1(G) C V(G) to be the set of vertices of G that are not cut vertices 
and additionally are not equal to 0. It is a simple fact that every connected multigraph 
on at least two vertices has at least two vertices that are not cut vertices. Thus T(G) 
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is in fact nonempty. The inductive step in the argument will be to consider G — i and 
its set of parking functions for some i 6 1(G). 

For any 1 < i < n let MPF(G; i) be the subset of a = (ai,.... a n ) £ MPF(G) for 
which ai is maximized. We claim that for i £ 1(G) the map /3 i—>• /3 + (deg G (i) — l)[z] 
is a bijection between MPF(G — i) and MPF(G;i) (when we continue to view 0 as 
the root of G — i). Why is this? First, note that an equivalent description of A(G) 
is as the set of acyclic orientations O of G for which there is a path from 0 to i for 
every i £ V(G). Here & path in O from s to t is a sequence (*2, *3), • • •, (ikAk+ 1) 

with i\ = s, 1 = t , and 1) £ O for all 1 < j < k. For i £ V(G) let A(G\i) be 

the set of those orientations O £ A(G) for which i is a sink. The alternate description 
of A(G) implies that for any 1 < i < n there is a bijection between A(G; i) and A(G — i) 
given by O i-A O D V(G — i) 2 . Next, note that A{G — i) = 0 if and only if i is a cut 
vertex (one way to see this is via Fact 12.41 we have MPF(G — i) 7^ 0 if and only 
if PF(G — i) 7^ 0 and, as mentioned earlier, PF(G — i) 7^ 0 if and only if G — i is 
connected). Thus A(G\i) is nonempty for all i £ 1(G). But as long as A(G;i) is 
nonempty, Fact 12.41 savs that the maximal value of a. % among (oq,... ,a n ) £ MPF(G) 
must be deg G (i) — 1. Moreover, in this case the elements of MPF(G) that achieve 
this maximal 07 are exactly the ones whose corresponding orientations (under the 
bijection of Fact 12.411 belong to A(G\i). That the map (3 i-A (3 + (deg G (i) — 1)[*] is a 
bijection between MPF(G — i) and MPF(Gr; i) for i £ 1(G) then follows from the fact 
that O i-A O n V(G — i) 2 is a bijection between A(G] i) and A(G — i) (together with 
Fact 12.411 . 

There are two important consequences of the construction in the previous paragraph. 
The first consequence is that 

(*) deg G (i) = deg G (j) for all i,j £ 1(G). 

Why does 0 hold? Because, as mentioned above, if i £ 1(G) then the maximal value 
of ai among (aq,..., a n ) £ MPF(G) is deg G (i) — 1; but since MPF(G) is invariant 
under © n , this maximal value must be the same for all 1 < i < n. The second 
consequence is that 

(a*) G — i is an a-cycle or G is an a-tree or G — i = K!f ,b for all i £ 1(G). 

Why does fen) hold? The invariance of MPF(G) under 6 n implies that MPF(G;i) 
is invariant under the subgroup of & n consisting of permutations that fix i (this sub¬ 
group is of course isomorphic to S n _i.) But then the bijection between MPF(G;i) 
and MPF(G — i) from the last paragraph implies that MPF(G — i) is invariant un¬ 
der & n - 1. By our inductive hypothesis we conclude (j*FJ). 

As we pointed out earlier, 1(G) is nonempty. Without loss of generality let us 
assume that n £ 1(G). By (j**j) we know that G — n is an a-tree or G — n is an a-cycle 
or G — n = Kn’ b . We will consider each of these cases (although not in that order). 
First we define a bit of terminology that will aid in our analysis of these cases. For a 
multigraph H and i.j £ V(H), we say that i is adjacent to j in H if cUfl-({i, j}) / 0. 
We say H contains a cycle of adjacent vertices if there is a sequence ii, i%,... ,ik, h+i 
of vertices for some k > 3 such that i\ = ik+i, ij and ij+i are adjacent for all 1 < 3 < k, 
and ij 7^ ij / for any 1 < j,j' < k. Note that if H is an a-tree then H does not contain 
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a cycle of adjacent vertices. Also note that if H is an a-cycle and #V(H) > 3 then 
every vertex of H is adjacent to exactly two vertices. And if H = then every 

vertex of H is adjacent to every other vertex. 

Case I: G — n is an a-cycle with a > 1 and n > 4. 

First suppose that n is adjacent to only one i € V(G ). We must have u ;©{*, n-}) = 2a 
because of © • Take any 1 < j < n with i ^ j and observe that G — j is not an a'-tree 
(because it has an edge of weight wg({*, n}) = 2a > a as well as an edge of weight a) or 
an a'-cycle or equal to I\n ,b (because it has a vertex adjacent to only one other vertex), 
contradicting Thus n must be adjacent to more than one vertex. Now suppose 

that n is adjacent to at least two other vertices. Then note that 1(G) = V(G) and so 
because of © we must have ujc({i,n}) = uJc({j,n}) for all 1 < i,j < n. Let 1 <i <n 
be such that i is adjacent to 0. We can see that G — * is not an a'-tree (because it 
contains a cycle of adjacent vertices) or an a'-cycle (because it has a vertex adjacent 
to at least three other vertices) or Kn ,b (because it has at least two vertices that are 
not adjacent), contradicting (j**j) . Therefore this case is actually impossible. 


Case II: G — n = Kn’ b with a, b > 1 and n > 3. 

This case includes when G — n is an a-cycle and n = 3. Suppose that n is adjacent to 
only one i £ V ( G ). Then u({i, n}) = a(n — 1) + b because of ©. Take any 1 < j < n 
with i 7^ j and observe that G — j is not an a'-tree (because it contains an edge of 
weight u({i, n }) = a(n — 1) +b > b as well as an edge of weight b ) or an a'-cycle or Kn ’ b 
(because it has a vertex adjacent to only one other vertex), contradicting (k*j) . Thus n 
must be adjacent to more than one vertex. Now suppose that n is adjacent to at least 
two other vertices. Then note that T(G) = V(G) and so because of © we must have 
UG({i,n}) = UG({j,n}) for all 1 < i,j < n. Let 1 < i < n. We see that G — i cannot 
be an a'-tree (because it contains a cycle of adjacent vertices) and G — i cannot be an 
a'-cycle with n > 4 (because any two of its vertices are adjacent). So from ©Fl) we 
conclude G — i = I\n ’ b . Moreover, by considering edges that belong to both G — i 
and G—n we must have a = a' and b = 6'. This implies in particular that wg({0, n}) = a 
and uic({j, n}) = b for all 1 < j < n. Thus G = j. 

Case III: G — n is an a-tree with a > 1. 

This case includes when G — n is an a-cycle and n = 2, as well as G — n = K? 2 ’ b . 
First suppose n is adjacent to only one vertex. Then A(G) has only one element and 
Fact 12.41 implies that the only way that MPF(G) is invariant under & n is if G is an 
a-tree. Thus, suppose from now on that n is adjacent to at least two vertices. If n = 2 
then we can see from © that G = K^’ b since 1(G) = {1,2}. So now suppose n = 3. 
First suppose that G — 3 is 


(t) 


» 

1 


a a 

—j— 

0 


2 


(Here we label an edge e £ E(G) by uc(e), with an edge absent if it has weight zero.) 
With G — 3 as in 0 we have 1(G) = {1,2,3}. Because of © it is impossible that 3 
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is adjacent to only 1 and 0, or to only 2 and 0. And because of (jvj) if 3 is adjacent to 
only 1 and 2 then G = Cf. Lastly, let us rule out the possibility that G looks like 



for some b,c,d> 1. By applying (j**j) to 1, we see that we must have b = a. But then G 
cannot satisfy Q because deg G (3) = a + c + d>a + c = deg G (l). So indeed it is not 
possible for 3 to be adjacent to 0, 1, and 2 when G — 3 is as in ([j]) . We have considered 
all possibilities for G — 3 as in ([]]). Next suppose that G — 3 is 

a a 

4-4 4 


(ft) 


0 1 2 


Because of (jvj) it is impossible that 3 is adjacent only to 1 and 0. And because of (j*j) 
if 3 is adjacent to only 0 and 2 then G = C%. Let us rule out the possibility that G 
looks like 



for some b,c,d> 1. By applying (TvaD to 2, we see that we must have b = a. But then G 
cannot satisfy (gj) because deg G (3) = a + c + d>a + d = deg G (2). So indeed it is not 
possible for 3 to be adjacent to 0, 1 , and 2 when as in ( Iff] ). Lastly, let us rule out the 
possibility that G looks like 



for some b, c > 1: if it did, then for all O € A(G) we have indeg G (l) = a < a + c 
but there is some O' E A with indeg G /(2) = a + c which by Fact 12.41 directly implies 
that MPF(G) is not invariant under & n . We have considered all possibilities for G — 3 
as in ©• Up to isomorphism, G — n can only look like ([[]) or fffj ) when n = 3. 

So finally let us suppose that n > 4 . Let T be the tree such that uic-n = cl ■ cut- 
Observe that X(G — n) is the subset of elements in {1 ,..., n — 1} that are maximal with 
respect to <t, so in particular Z(G — n) is nonempty. Note also that X{G — n) C X{G) 
because n is adjacent to at least two vertices. First suppose that n is adjacent to 
three or more vertices. Then for any i € X[G — n) we will get that G — i is not Kn ,b 
(because it contains two vertices that are not adjacent) and is not an a'-tree (because 
it contains a cycle of adjacent vertices). Moreover, by our analysis of Case I we also 
know that G — i cannot be an a'-cycle for any i E 1(G). This contradicts tell . So it 
is not possible that n is adjacent to three or more vertices. Next suppose that there 
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is i E X{G — n) which is not adjacent to n. Then we will again have that G — i is 
not Kn ,b (because it contains two vertices that are not adjacent) and is not an octree 
(because it contains a cycle of adjacent vertices) and is not an a'-cycle (by appealing 
to our analysis of Case I), contradicting (j**j) . We conclude that n is adjacent to 
exactly two vertices; more specifically, n is adjacent to all of the vertices in X{G — n ) 
(which contains either one or two vertices) and at most one other vertex (in the case 
where j^X[G — n) = 1). Therefore, because I(G — n) is the set of maximal elements 
of T with respect to <T, and also by applying (j**j) to each i E X(G) to check that G — i 
is an a-tree, we see that G looks like 



(The dotted lines indicate paths of edges of weight a.) Suppose vertex j above is 
not 0: then for all O E A(G) we have indeg^j) = a < 2a but there is O' E A 
with indeg£>/(?r) = 2 a which by Fact 12.41 directly implies that MPF(G) is not invariant 
under & n . Thus we must have that j = 0 and that G is an a-cycle. 

To finish the proof, it is straightforward to verify that in all these cases where PF(G) 
is invariant under & n it is equal to PF(x) for the claimed x = (aq,... ,x n ) E N”. For 
this purpose it is useful to observe that a E N n is an x-parking function if and only 
if a < a! for some maximal x-parking function a!. Thus, thanks to Fact 12.31 we 
need only check that the maximal x-parking functions are the same as the maximal 
G-parking functions for the appropriate G and x. The maximal G-parking are easy to 
understand because of Fact 12.41 while the maximal x-parking functions are merely the 
permutations of (aq — 1, x\ + aq — 1 ,..., x\ + X 2 + ■ ■ ■ + x n — 1). □ 

Remark 2.6. There are simple product formulas for ^PF(x) for all of the x appearing 
in Theorem 12.51 Namely, 

n— 1 

. #PF((a,67C _ 0)) = a"; 

n —2 

. #PF((a,oXX0,a)) = (n + lK; 

n— 1 

• #PF {{a,wC~~b) = a{a + nb) 11 - 1 . 

The first two formulas are trivial; for the third see [43l Equation (7)] or [1471 Theorem 1]. 

3. The multigraph DFS-burning algorithm 

In this section we extend the work of [33] to multigraphs. Thus for a multigraph G we 
relate the reversed sum enumerator of G-parking functions to the ^-number (generalized 
inversion number) enumerator of spanning trees of G. We need to define these k- 
numbers. Let G be a fixed connected multigraph with V(G) := {0,1,...,n} and 
root 0. Following Gessel [Tfi] . for a spanning tree T E SPT(G) and a total order -< 
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on V (G) we define 


k {G,T, -<) := ^ w G ({par T (i),j}). 

i,j£V(G), 
i<rj, 3<i 


Theorem 3.1. We have 


k(G,T,~<) . 


n Me)], 

ye&E(T) 


^ ^rsum(a) _ 

aePF(G) TeSPT(G) 

where -< is any total order on V(G). 

Theorem 13.II will follow immediately from the existence of a bijection 


</? . PF(G) —>• < (T,£): 


Tg SPT(n + l), £: E(T) -A N, 

£(e) G {0,1,... ,wc(e) — 1} for all e G E(T) 


satisfying rsum(ai) = k(G, T, -<) + Yle£E(T) ^( e ) when (T, (a). We establish the 

existence of such a bijection in Theorem 13.31 below. This bijection is a straightfor¬ 
ward extension of the DFS-burning algorithm of m to multigraphs. 

Let us recap the “generic Dhar’s burning algorithm.” Dhar’s burning algorithm 
takes an input a G-parking function and outputs a spanning tree of G. So suppose we 
are given some a = (au, ... ,a n ) G PF(G). To initialize, we put ctj chips on vertex i 
for all 1 < i < n. The process is called the “burning algorithm” because we imagine 
that a fire is spreading throughout the graph: the chips can be seen as impediments 
that the fire has to burn through. At each step of the burning algorithm we have 
some collection of vertices that are burning together with a rooted tree that spans the 
burning vertices. Initially only the root 0 is burning. At each step of the algorithm 
we choose an edge connecting a burning vertex to one that is not yet burning. If the 
non-burning vertex has a positive number of chips on it, we remove one chip from that 
vertex and we burn the edge we selected, removing it from consideration at later steps 
of the algorithm. If the non-burning vertex has no chips on it, that vertex starts to burn 
and we include the selected edge in the tree we are building up. Once all the vertices 
are burning we terminate and output the resulting spanning tree T. Importantly, the 
burning algorithm keeps track of the sum of its input: the sum of a is precisely the 
number of edges that were burnt. The reversed sum is the number of edges that were 
not burnt and are not in T; we call these edges the surviving edges. 

To turn the “generic” burning algorithm into an actual, deterministic process we 
need to specify a procedure for choosing an edge between a burning and non-burning 
vertex at each step: let us call this choice the choice of which edge to burn along. 
In the Cori-Le Borgne |T0] variant of the burning algorithm we have some fixed total 
order on the edges and always choose to burn along the maximum edge according to 
this order. The surviving edges are precisely the edges that are externally active in 
the output tree T; consequently, the Cori-Le Borgne algorithm relates reversed sums 
of G-parking functions to external activities of spanning trees of G. In the DFS-burning 
algorithm of Perkinson-Yang-Yu [35] we have some fixed total order -< on the vertices 
and choose to burn along the edge connecting the burning vertices to the maximum 
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non-burning vertex according to -<, but in a depth-first search fashion. In the case of a 
simple graph G (i.e., a graph with no multiple edges), the surviving edges are precisely 
those of the form {par T (i),j} with i <t j but j -< i (where T is the output tree); 
consequently, the DFS-burning algorithm relates reversed sums of G-parking functions 
to K-numbers of spanning trees of G. Incidentally, there are many other possible graph 
search procedures that one could use to specify which edge to burn along. Some of 
these, such as breadth-first search and the neighbor-first search of Gessel and Sagan m 
were explored by Kostic and Yan [25). A large family of graph search procedures were 
also investigated by Chebikin and Pylyavskyy in [9]. Each choice of graph search 
procedure relates reversed sum to some new spanning tree statistic. 

The one wrinkle that can occur when applying the DFS-burning algorithm to multi¬ 
graphs is that some of the surviving edges can be parallel to edges of the output tree T, 
in which case these edges do not correspond to inversions. But all we have to do to 
correct for this is to give the output tree T an edge-labeling function £: E(T ) —>• N that 
records how many surviving edges were parallel to each e £ E(T). Algorithm [T| gives 
pseudocode for the resulting multigraph DFS-burning algorithm. Algorithm [2] gives 
pseudocode for the inverse of this algorithm: the inverse algorithm uses a very similar 
depth-first search burning procedure but now instead of removing chips, it adds chips 
to vertices as edges are burnt. 


Algorithm 1 Multigraph DFS-burning algorithm. 

ALGORITHM (WITH RESPECT TO MULTIGRAPH G AND TOTAL ORDER -< ON V(G)) 
Input: a = (an,... ,a„) £ N n 

1: burning_vertices := {0} 

2: burnt_edges := 0 

3: T ■= multigraph with V(T) := {0} and E(T) := 0 
4: £ := unique map E(T ) -A N 
5: head := 0 

6: foreach j £ V(G) — burning_vertices in order from max to min according to -< do 
7: foreach e = {head, j} £ E(G) — burnt_edges do 

8: if aj > 1 then 

9: burnt_edges := burnt_edges U {e} 

10: aj := otj — 1 

11: else 

12: burning_vertices := burning_vertices U {j} 

13: V(T) := V(T) U {j} 

14: E(T) := E(T) U {e} 

15: £{e) := #{/ £ E{G) — (burnt_edges U {e}) : / = {head, j}} 

16: head := j 

17: break out of current for loops and go to line [G] 

18: if head 0 then 
19: head := par T (head) 

20: go to line [6] 

Output: (T,£) 
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Algorithm 2 Multigraph DFS-burning inverse algorithm. 

ALGORITHM (WITH RESPECT TO MULTIGRAPH G AND TOTAL ORDER -< ON V(G)) 

Input: (T, t) with T G SPT(G) and £: E(T) -» N 
1: burning_vertices := {0} 

2: burnt_edges := 0 

3: a = (ai, ..., a n ) := (0,0,..., 0) 

4: head := 0 

5: foreach j G V(G) — burning_vertices in order from max to min according to -< do 
6: foreach e = (head, j} G E(G) — burnt_edges do 

7: if e G E(T) and £(e) = #{/ G E(G) — (burnt_edges U {e}) : / = {head, j}} then 

8: burning_vertices := burning_vertices U {j} 

9: head := j 

10: break out of current for loops and go to line [5] 

11: else 

12: burnt_edges := burnt_edges U{e} 

13: otj '. = Oij 1 

14: if head ^ 0 then 
15: head := par T (head) 

16: go to line [5] 

Output: a 


Example 3.2. Let G be the multigraph below, where for clarity we have labeled each 
edge e G E(G) by wc(e) and each vertex i G V(G) by vf 



Let a := (1,2,0,0) G PF(G) and let -< be given by 1 -< 2 -< 3 -< 4. Figure [2] depicts a 
run on the multigraph DFS-burning algorithm on a with respect to the vertex order 
The result of the DFS-burning algorithm is that (T,£) := (p^(a) is 



where now we label each e G E(T) by £(e). We can verify that rsum(a) is 
4 = wg({0, 2}) + cjg({0, 3}) + cjg({4, 1}) + cjg({0, 1}) + 2 = k(G,T, -<) + ^(e). 

eG E(T) 

Observe also that 4 is the number of surviving edges, that is, the number of solid edges 
not in T in the lower-left of Figure [21 
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Figure 2. Example 13.21 a run of the multigraph DFS-burning algo¬ 
rithm. The algorithm starts in the upper-left and ends in the lower-left. 
Edges of the graph are drawn with multiplicity for clarity. At each step, 
the burning vertices are marked by an asterisk and the non-burning 
vertices are marked by the number of chips remaining on them. The 
head of the depth-first search has a larger asterisk. The edges in the 
tree spanning the burning vertices are bold and directed away from the 
root; their labels are displayed beside them. The burnt edges are dashed. 


Theorem 3.3. Let (T. £) be the output of an application of Algorithm, [7] to a £ N n 1 
with respect to any total order -< on V(G). IfV(T) = V(G) then a € PF(G'). And 
if V(T) / V(G), then the set U := V{G) — V(T) certifies that a is not a G-parking 
function in the sense that for all j £ U we have aj > dy(j). 

For a fixed total order -< on V(G), associating to each a £ PF(G) the edge-labeled 
spanning tree (T. t) produced by Algorithm^ on input a yields a bijection 


: PF(G) —» < (T,£): 


re SPT(G), l\ E(T) -A N, 

£(e) £ {0,1,... , wc(e) — 1} for all e £ E(T) 

This bijection satisfies rsum(o;) = k(G, T, -<) + ^eeE(T) ^( e ) w ^ en {T,£) = The 

inverse of the bijection is given by Algorithm [H 


Proof. The proof of the correctness for this algorithm is so similar to that of the original 
algorithm in [35] that we will not repeat it here. □ 

n —1 

Remark 3.4. When x = (a, b, b,... , b) we have PF(x) = PF(K by Theorem 12.51 
Thus Theorem 13.11 recovers a result of Yan [41] expressing the reversed sum enumerator 
for these x-parking functions in terms of inversions in “rooted 6 -forests.” 
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4. The vector DFS-burning algorithm 


In this section we prove our main result, Theorem 11.81 Fix x = (xi,... ,x n ) £ N n . 
Assume xi > 1 as otherwise Theorem 11.81 holds trivially since PF(x) = 0. Define a 
connected multigraph G x associated to x: its vertex set is V(G x ) := {0,1 ,..., n} and 
its edge-weight function is given by wg x ({*>J'}) := x 'min(?jj+i for all 0 < i,j < n. As 
always, 0 is the root of G x . Note that #E(G X ) — n = (Xa= i( n + 1 — i)xi) — n so 
the notions of reversed sum for G x -parking functions and x-parking functions agree. 
Theorem Ol will follow immediately from the existence of a bijection 


ip : PF(x) —> 




T £ RPT(n + 1), -<£ AVO(T), £: E(T) -A N,1 
£{e) £ {0,1,... ,WG x (e) — 1} for all e £ E(T) J 


satisfying rsum(a) = k(G x ,T, -<) + Yl e eE(T) ^( e ) when (T,£, -<) = ip (a). We establish 
the existence of such a bijection in Theorem 14.21 below. This bijection is a far more 
substantial variation of the DFS-burning algorithm than the extension to multigraphs 
in the last section. We want to “symmetrize” the DFS-burning algorithm. The trick 
is that we will build up our graph G x as we burn through it. We start with a graph 
on {0,1,..., n} which has no edges. Then 0 starts to burn. When a vertex j becomes 
the ith vertex to start burning we add Xi edges between j and all vertices that are not 
yet burning. An auxiliary permutation a in the algorithm records the order in which 
vertices started to burn. The permutation a determines the admissible vertex order -< 
we return together with the rooted plane tree T and its edge-labeling function £. When 
we terminate the algorithm the graph we built up will always be a copy of G x once 
we relabel the vertices by * e> u _1 (i). The point is that this vector DFS-burning 
algorithm simulates a run of the multigraph DFS-burning algorithm on this relabeled 
multigraph G x with respect to the order -< we obtain as output. Therefore, when 
analyzing this new procedure we can invoke all the nice results about the behavior of 
the multigraph DFS-burning algorithm detailed in Theorem 13.31 that it respects the 
reversed sum statistic, that it is invertible, et cetera. Algorithm [ 3 ] gives pseudocode 
for the vector DFS-burning algorithm. Algorithm [4] gives pseudocode for the inverse of 
this algorithm, which just applies the multigraph DFS-burning inverse algorithm. 


Example 4.1. Set x := (1,3,1) and take a := (2,0,2) £ PF(x). Figure [3] depicts a 
run of the vector DFS-burning algorithm on a. As in Example 13.21 for clarity we label 
vertex i £ V(G) by Vj. But note that the output tree T will have its vertices relabeled 
by i <L~ 1 (*) with respect to the depicted Vi labeling in Figure [3J Specifically, the 
result of the DFS-burning algorithm is that (T, £, -<) := ip(a) is 


rc 0 



T 


where we label each i £ V(T) by Wi, each e £ E{T) by £(e), and 3 -< 1 -< 2. We can 
verify that rsum(a) = 3 = wg x ({0, 3}) + 2 = n(G, T, -<) + ^eeE(T) ^( e )- 
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Algorithm 3 Vector DFS-burning algorithm. 

ALGORITHM (WITH RESPECT TO X = {x\,...,X n ) £ N”) 

Input: a = (an,... ,a„) £ N n 

1: G := multigraph with V (G) := {0,1,..., n} and E{G) := 0 
2: burningjvertices := {0} 

3: burnt_edges := 0 

4: counter := 0 

5: a := unique map {0} —> {0} 

6: foreach k £ V(G) — burning_vertices do: 

7: Add x\ copies of the edge {0, k} to E(G) 

8: T := multigraph with V{T) := {0} and E(T) := 0 
9: l := unique map E(T) -A N 
10: head := 0 

11: foreach j £ V(G) — burning_vertices in order from max to min j do 
12: foreach e = {head, j} £ E(G) — burnt_edges do 

13: if Oij > 1 then 

14: burnt_edges := burnt_edges U {e} 

15: ctj := ctj — 1 

16: else 

17: burning_vertices := burning_vertices U {{} 

18: counter := counter + 1 

19: er(counter) := j 

20: foreach k £ V(G) — burningjvertices do: 

21: Add a.’counter+i copies of the edge {j, A;} to E(G) 

22: V(T) := V(T) U {j} 

23: E(T) := E(T) U {head, j} 

24: f({head, j}) := #{/ £ E(G) - (burnt.edges U {e}): / = {head,.;}} 

25: head := j 

26: break out of current for loops and go to line fill 

27: if head ^ 0 then 

28: head := par T (head) 

29: go to line [IT] 

30: relabel the vertices of T by m ct _1 (j) for all i £ V(T ) 

31: := the unique total order on V(T) with i -< j if and only if a (i) < a U) 

Output: (T, £, -<) 


Algorithm 4 Vector DFS-burning inverse algorithm. 

ALGORITHM (WITH RESPECT TO X = (xi,. . . , X n ) £ N”) 

Input: (T, e, -<) with T £ RPT(n + 1), £: E{T) -A N, and ^!£ AVO(T) 

1: /3 = (/3i,. .., j 3 n ) := result of applying Algorithm [I] to (T, €} with respect to G x and -< 
2: a \= the unique permutation in 6 n with <r(l) -< <r(2) -<•••-< cr(n) 

^ • (/^ct(1) ? Pcr(2) ) • * * ) Pcr(n) ) 

Output: a 
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V 3 v 2 ,a(l) 


V 3 V 2 ,<t(1) 





v 0 ,cr(0) 


V 3 ,a( 2 ) v 2 ,<r(l) v 3 ,a(2) v 2 ,cr(l) v 3 ,a(2) v 2 ,a(l) v 3 ,a(2) v 2 ,a(l) 





Figure 3. Example 14.11 a run of the vector DFS-burning algorithm. 
The algorithm starts in the upper-left and ends in the lower-left. 


Theorem 4.2. Associating to each x-parking function a £ PF(x) the triple (T,£,~<) 
produced by Algorithm [5j on input a defines a bijection 


if: PF(x) —> 


(TJ, ~<) 


T £ RPT(n + 1), -<£ AVO(T), £: E(T) -A N,1 
£(e) 6 {0,1,... ,cu Gx (e) - 1} for all e £ E(T) } 


This bijection satisfies rsum(a) = k(G x ,T, -<) + YleeEtT) ^( e ) when (T,£,~<) = if (a). 
The inverse of the bijection is given by Algorithm [7} 


Proof. Let a £ PF(x) and let (T,£, -<) be the output of the vector DFS-burning algo¬ 
rithm on input a. First let us show that indeed T £ RPT(n +1) and A£ AVO(T). It is 
clear that T is a tree. We need to show that in fact T has n +1 vertices. Suppose to the 
contrary that i := ffV(T) < n + 1. Let G,a be as in Algorithm [3] when the algorithm 
terminates. Set U := V(G) — cr(V(T)). All edges in G between U and a(V(T)) were 
burnt without any vertex in U starting to burn, which means that for each j €. U we 
have ctj > deg jj(j) = J2j =o x i• Thus there are at least n — i entries of a greater than or 
equal to J2j=o X T means that a is not an x-parking function, a contradiction. 

So indeed V(T) = {0,1,..., n}. We have T £ RPT(n + 1) because the algorithm adds 
vertices to T in a depth-first search fashion; the vertices of T have been relabeled in 
such a way that vertex i was the (i + l)st vertex added to T. We have -<£ AVO(T) 
because the algorithm prefers to visit vertices j with a greater label first. 

Continue to fix some a € PF(x) with (T,£, -<) := if (a). Let a £ & n be the unique 
permutation with <j(l) -< a(2) -<•••-< a(n). Set fi := ot a -i( 2 )> ■ ■ ■ j a <r-i(n))- 

The crucial observation is the following: a run of the vector DFS-burning algorithm with 
input a simulates a run of the multigraph DFS-burning algorithm (Algorithm [T|) with 
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input f3 with respect to multigraph G x and vertex order -<. Thus, because the algorithm 
terminates with V(T) = V(G X ), Theorem 13.31 implies that j3 G PF(G X ). Moreover, we 
have that (T ,£) = <p^((3). One consequence of (T ,£) = f^(/3) is that £: E(T ) -A N 
satisfies £(e) G {0,1,... ,wc x (e) — 1} for all e G E(T). We conclude that ^ is well- 
defined. Another consequence of (T ,£) = (p^(/3), again appealing to Theorem 13.31 
is 

rsum(a) = rsum(/3) = k(G x , T, A) + E 

eSE(T) 

Furthermore, Theorem [373] implies that if is injective with inverse given by Algorithm^! 

To finish the proof we need to show that if is surjective, or equivalently that when we 
apply the inverse algorithm to a valid triple (T, £, -<) we obtain an x-parking function. 
So let T G RPT(n+l), AG AVO(T), and C E(T) -A N with £(e) G {0,1,... , w Gx (e)-l} 
for all e G E(T). Let a G N n be the output of the inverse algorithm (Algorithm [5j) 
on input (T, £, -<). Let (3 be as in Algorithm |4j Because T has a depth-first search 
labeling, and because A is an admissible vertex ordering, the order that vertices of G x 
start to burn as we carry out the inverse algorithm is precisely 0,1,2, ...,n. Thus 

for any 1 < i < n, at most deg |( x +1 n j(i) — 1 = — 1 chips were added 

to i as we carried out the inverse algorithm. This means that f3i < x i ) ~ 1 

for all 1 < i < n, which implies that (3 G PF(x). But a is a permutation of (3 and 
therefore a G PF(x) as well. □ 

As mentioned, Theorem 11.81 follows immediately from Theorem 14.21 Now let us see 
why Corollary 11.111 follows from Theorem 11.81 

Proof of Corolla,rv \l.ll\ It suffices to show ffAVO (T) = n! 11^=1 (outdeg^^ — l )!) -1 for 
all T G RPT(?r + 1). To see this, let T G RPT(n + 1). A total order A on V(T) is 
admissible if and only if for each i = 0 , 1 ,,n the children of i appear in decreasing 
order of label in A. Clearly (outdeg r (i )!) _1 °f the n! total orders on V(T ) have the 
children of i in the appropriate order. Moreover, these conditions are all independent 
since each vertex has only a single parent. Lastly, note that n cannot have a child in T 
so we can omit the term (outdeg T (n )!) -1 from the product. □ 

Remark 4.3. Let us show how Corollary II. Ill is equivalent to Theorem [L5] in a simple 
way. The point is that T i-a y(T) := (outdeg T (0), outdeg r (l),..., outdeg T (n — 1)) 
defines a bijection between RPT(n + 1) and T(n). Why is this? Let T G RPT(?r + 1). 
First of all, 7 (T) is indeed in T(n) since 

• Xa=i t(-O* > j for 1 < j < n because par T (i) G {0,1,..., j — 1} for 1 < i < j\ 

• Sr=i 7(^)i ^ n because only {1,..., n} are children in T. 

Moreover, it is not hard to see that for any i G V{T ) the set of descendants of i in T 
is precisely {i + 1 , i + 2 ,..., i + m} where 


m := ma x{k G {0,1,... ,n-i }: ^{T) i+1 +^{T)i +2 3 -by {T) i+ j > j for all 1 < j < k}. 
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But knowledge of the descendants of i for each i E V(T) determines the tree T. So 
the map T H y 7 (T) is injective. That T 1 —> 7 (T) is bijective then follows from the 
well-known facts, mentioned in Section [I] that #RPT(n + 1) = C n = 


5. Open questions and future directions 
The following are some possible threads of future research. 

(1) It is often worthwhile to substitute q := 1,0, —1 in any statistical enumerator for a 
combinatorial set. Specializing q := 1 in the reversed sum enumerator for PF(x) of 
course yields #PF(x). Specializing q := 0 counts the number of maximal x-parking 
functions: these are just permutations of the partition 


A x := (x n + x n -i -I-f- x\ — 1,..., x 2 + x\ — 1, x\ — 1) 

and so are easily enumerated. Chebikin and Postnikov [8] (see also Pak and Post- 
nikov [33]) showed that the absolute value of the specialization q := —1 is 0 if A x is 
even and is /3 n (S x ) if A x is odd. Here /3 n (S ) for S C { 1 , 2, ..., n— 1} is the number of 
permtutations in & n with descent set S, and S x := {i: Af +1 is odd, 1 < i < n — 1}. 
Via Theorem 11.81 the q := 1 specialization yields the expression for ^PF(x) in 
terms of trees stated in Corollary 11.111 It would be interesting to interpret both 
the q := 0 and q := —1 specializations in terms of trees via Theorem ll.81 The q := 0 
specialization should lead to an extension of the notion of increasing tree (see 
e.g. HD Proposition 1.5.5]). Meanwhile, one could hope to find a involution on 
rooted plane trees with admissible vertex orders that inverts A-number parity and 
recovers, by consideration of fixed points, the aforementioned formula for q := — 1. 

(2) A weakly increasing x -parking function is a = («i,...,a n ) E PF(x) that satis¬ 
fies 07 < a .2 < ■ ■ ■ < a n . Denote the set of weakly increasing x-parking functions 
by PF^x). It is not hard to see that weakly increasing x-parking functions are 
in bijection with partitions ji = with /u C A x (here n C A x means 

that Hi < A x for all 1 < i < n). Moreover, this bijection yields the identity 


E 


^rsum(a) 


aePF/'fx) 


q \^\-\u\ 

/iCA x 


where for a partition v = ( 17 ,..., u n ) the size of u is \v\ := D- It would be 

interesting to give an expression for the reversed sum enumerator for PF^'(x) in 
terms of trees and inversions. The DFS-burning algorithm does not respect weak 
increasingness in any obvious way. Specializing q := 1 of course yields ^PFP'(x) 
or equivalently the number of partitions h contained in a given partition A x . This 
number is given by a determinant that was essentially known to MacMahon 130 
p. 243]; see the discussion after Theorem 12 in (43j . 

(3) Although Theorem 12.51 classifies the graphs whose parking functions are also vector 
parking functions, there can still be exceptional equalities between the numbers of 
parking functions for other families of graphs and vectors. Here is an example of 
such an exceptional equality. Let I \^ m denote the complete bipartite graph with 
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edge-weights a; e.g. V(K^ m ) := { 0 , 1 ,, 2 m — 1 } and 

., ,. la if 7 — i = 1 mod 2 

uj K a ({z,j}) := < 

ra ’ mU S> |0 otherwise. 


Set y (a, 2, m) := (a, 0, a, 0,... , a) € N 2m D Then for all a > 1, m > 2, 

#PF(A- ;m ) = a 2m - l m 2m ~ 2 = #PF(y(a, 2,m)). 

The map a o 2a bijects between PF(y(l,2, m)) and the subset of PF(2m — 1) 
for which all entries are even, so we might call PF(y(l, 2, to)) the set of “even 
parking functions.” We can slightly extend this example by considering directed 
graphs. There is a notion of G-parking function for G a directed graph (see HE [HD- 
Let G(a, b , m) be the directed multigraph with V ( G(a , b, to)) := {0,1,..., bm — 1} 
and with directed edge-weight function 


^G(a,b,m) (jt f) 


a if j — i = 1 mod b 
0 otherwise. 


Note that G(a, 2 , to) = m and G(a, 1, to) = A'S“ are actually undirected graphs. 
Set 


y(a, 6 , to) 



6-1 6-1 6-2 

✓s. 

, 0 ,... , 0 , a, 0 , 0 ,... , 0 ,..., a, 0 , 0 ,... , 0 ) € N 


bm— 1 


Then for all a, b > 1 and to > 2, 

#PF (G(a,b,m)) = a brn ~ 1 m bm ~ 2 = #PF(y(a, 6 , to)). 

Apparently PF(G(a, 6 , m)) 7 ^ PF(y) for all choices of a > 1, 6 > 2, and m > 3. 
It is likely possible to extend the classification of graphs whose sets of parking 
functions are invariant under & n in Theorem 12.51 to include directed graphs but 
such a classification is beyond the scope of the present paper. At any rate, it 
would be very interesting if the equality #PF (G(a,b,m)) = #PF(y(a, 6 , to)) were 
more than just a numerical coincidence: is there some reasonable bijection between 
the two sets? Even in the case a := 1 and b := 2, we are not aware of a bijection 
between even parking functions and spanning trees of the complete bipartite graph. 
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